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Two-pion correlation functions in Au+Au collisions at -^/sjvjv = 130 GeV have been measured 
by the STAR (Solenoidal Tracker at RHIC) detector. The source size extracted by fitting the corre- 
lations grows with event multiplicity and decreases with transverse momentum. Anomalously large 
sizes or emission durations, which have been suggested as signals of quark-gluon plasma formation 
and rehadronization, are not observed. The HBT parameters display a weak energy dependence 
over a broad range in y / sjvjv- 



Two-particle intensity interferometry techniques 
(HBT) have been used extensively to probe the space- 
time structure of heavy ion collisions . At midrapidity 
(y = 0) and low transverse momentum (j>t), two-pion 
correlation functions reflect the space-time geometry of 
the emitting source, while dynamical information (e.g. 
collective flow) is contained in the momentum depen- 
dence of the apparent source size J2|-||] . 

In this Letter we study two-pion correlation functions 
in Au+Au collisions at nucleon-nucleon center-of-mass 
energy ^/snn — 130 GeV produced by the RHIC facility 
at Brookhaven National Laboratory, and compare our re- 
sults to similar studies at lower energy. Such systematic 
comparisons may reveal the onset of new phenomena, as 
^Jsmn increases. We are particularly interested in pos- 
sible indications that a quark gluon plasma (QGP) has 
been formed in the collision. Several authors flH-pf have 
proposed HBT studies to probe for a significant increase 
in the pion emission timescale associated with QGP for- 
mation. 

Experimentally, the two-particle correlation function is 
obtained from the ratio Ca(q) = A(c\)/B(q) (normalized 
to unity at large q), where A(q) is the measured two-pion 
distribution of pair momentum difference q = P2 Pi , 
and B(q) is the mixed background distribution 0, cal- 
culated in the same way using pairs of particles taken 
from different events. 

The STAR time projection chamber (TPC) 0,|ll|] was 
used to record charged particle tracks as they crossed up 
to 45 padrows of the detector. For P-10 gas in the 0.25 T 
solenoidal magnetic field, transverse and longitudinal dif- 
fusion were both about 350/xm times the square root of 



the drift length. The TPC was read out with 138,000 
waveform digitizer channels |12[ | ; the signal shaping width 
was comparable to the diffusion. For a typical 1.5-meter 
drift, two-track resolution was limited by the resulting 6 
mm signal spread, in addition to geometric spreading due 
to tracks traversing the padrows at finite crossing angles. 
This analysis used tracks with 0.125 < pr < 0.45 GeV/c 
and \y\ < 0.5, for which the reconstruction and it identi- 
fication efficiencies are high. We present results for three 
p T ranges: 0.125 - 0.225 GeV/c, 0.225 - 0.325 GeV/c, 
and 0.325 - 0.45 GeV/c. 

We select events with a collision position that is within 
75 cm of the mid-plane of the 4-m long TPC. An impor- 
tant feature of two-particle correlation functions is that 
single-particle phase space and acceptance effects cancel 
to first order. To preserve this feature in the collider 
environment (in which the collision vertex- and hence 
the single-particle acceptance- varies event-to-event), we 
only mix events which have a primary vertex longitudi- 
nal location within 15 cm of each other; however, mixing 
events over the full ±75 cm range introduces no signif- 
icant distortion in the correlation fit parameters, in the 
present analysis. 

Starting with a minimum-bias trigger p3[ , we char- 
acterize the centrality of collisions off-line according to 
the measured multiplicity of negatively-charged parti- 
cles with pseudorapidity \-q\ < 0.5. After accounting 
for event reconstruction inefficiency for low-multiplicity 
events, we estimate that the minimum-bias distribution 
contains ^90% of the hadronic Au+Au cross section p3[. 
We analyze ~ 10 5 events in each of three bins; bin 3 con- 
tains the 12% most central of the measured collisions, bin 



2 the next 20%, and bin 1 the next 40%. 

Particle identification was achieved by correlating the 
magnetic rigidity of a track with its specific ionization 
(dE/dx) in the gas of the TPC. In the momentum region 
of interest, n, K, and p/p are well-separated; contamina- 
tion of the pion sample by electrons is the primary con- 
cern. Based on simulations and on extrapolation from 
regions of clear e~ /n~ separation, we estimate that elec- 
trons comprise about 10% (4%) of the selected tracks in 
our lowest (highest) px bin. Tracks used in the correla- 
tions are required to project back to the primary interac- 
tion vertex within 3 cm, thereby selecting primarily pions 
emitted directly from the collision. Contamination from 
non-pions and pions originating from long-lived decays 
(e.g. A, oj) reduces the strength of the correlation, char- 
acterized by the parameter A, while leaving the extracted 
source radii unchanged (l4|,|l5| • 

Although track-splitting (incorrect reconstruction of a 
single particle as two particles) in the STAR TPC is a 
small effect in general (< 1%), it may have a strong im- 
pact on correlation studies at low \q\. False pairs gener- 
ated by track-splitting are removed by a topological cut 
which requires valid and distinct signals of both tracks on 
several TPC padrows. For the present analysis, residual 
effects of track-splitting are negligible. 

To eliminate the effect of track-merging (in which two 
tracks with similar momenta are reconstructed as a sin- 
gle track), we require the tracks in a pair to be well sep- 
arated (> 2.5 cm at the radius of the TPC inner field 
cage, 50 cm). In applying this cut also to the mixed- 
pair background, the effect of event-to-event variation in 
primary vertex position is taken into account. The anti- 
merging cut removes significant detector effects, but also 
discriminates against low-|g| pairs, which carry the cor- 
relation signal. This leads to an artificial reduction in 
the HBT parameters, which we estimate, based on sim- 
ulations flq| , to be 3-6% for the radii and 6-14% for A, 
depending on px- We correct for this effect in the results 
presented. 

We apply to each background pair a Coulomb correc- 
tion M corresponding to a spherical Gaussian source of 
5 fm radius; this correction is identical to that used by 
several other groups |14] , p0 21 j. In principle, this pro- 
cedure over-corrects the correlation function for realis- 
tic sources (e.g. in which some pions originate far from 
the core source). We find that the HBT radii change 
smoothly by ~ 10% as the strength of the Coulomb 
correction is varied from the standard correction to no 
correction. Here, we restrict ourselves to the standard 
Coulomb correction, allowing for a uniform extension of 
existing interferometry systematics from lower energy. 

One-dimensional correlation functions con- 
structed in the invariant quantity Qinv = 
\/(pi — P2) 2 — (Ei — E2) 2 are fit to the functional form 



C(Q mv ) = 1 + Aexp(-QL,iO 



(1) 



The 1-D correlation function (and fit) corresponding to 
our highest multiplicity class is shown in Fig. la for low- 
Px 7r~ . The 1-dimensional fit fails in lowest Qi nv bins, 
as has been observed in previous measurements p3,E3|. 
Hence, the extracted radius Ri nv ~ 6.3 fm is only a rough 
indication of the space-time extent of the source. 

A more detailed characterization of the emitting source 
is obtained through multidimensional correlation func- 
tions. We decompose the momentum difference q accord- 
ing to the Pratt-Bertsch [||,[7| "out-side-long" (indicated 
by o, s, and I subscripts) parameterization. Here, qi is 
parallel to the beam, q s is perpendicular to the beam 
and to the total momentum of the pair, and q is per- 
pendicular to qi and q s . Data are analyzed in the longi- 
tudinally co-moving system (LCMS) frame, in which the 
longitudinal component of the pair momentum vanishes. 
Figs, lb-d show one-dimensional projections of the cor- 
relation function C(q , q s ,qi) onto the q , q s , and qi axes, 
for 7r~ from the most central collisions. 

The 3-dimensional correlation functions are fit with the 
standard Gaussian form 

C(q ,q s , qi) = l + \exp(-R 2 <£ _ R 2 q 2 _ ^2^ (2) 

where Ri is the homogeneity length in the i direction [n| . 
Projections of the fit to the central collisions, weighted 
according to the mixed-pair background, are shown as 
curves in Fig. pi 

Systematic errors on the HBT fit parameters have two 
sources of roughly equal magnitude: (1) the uncertainty 
on the correction for the anti- merging cut, estimated 
equal to the correction itself, and (2) the uncertainty 
associated with the Coulomb correction, ~0.1-0.2 fm in 
the radii, determined by varying the Coulomb radius by 
±1 fm. 

The effect of the single-particle momentum resolution 
(Sp/p ~ 2% in the TPC for the particles under study) is 
known to induce systematic underestimates of HBT pa- 
rameters. For our lowest px bin, Sq « 4.5 MeV/c, and 
Sq s « 5qi sa 3 MeV/c. While Sq s and Sqi change little, in 
our highest px bin, Sq ss 11 MeV/c. Using an iterative 
procedure similar to that used at lower energies |14] , p2[ , 
we have corrected our correlation functions for finite res- 
olution effects. The correction had no effect on the HBT 
parameters for our lowest px cut; for the highest px cut, 
it resulted in a 5% increase in A and a 8% increase in R , 
while R s and Ri were unaffected. 

The multiplicity dependence of the source parameters 
is presented in the left panels of Fig. 0, for low-p^, midra- 
pidity 7r~7r~ and tt + tt + correlations. The parameters for 
positive and negative pions are similar; the A parame- 
ter is constant, and all three radii increase monotonically 
with multiplicity. The increase of the transverse radii 
R and R a with multiplicity is interpreted as a geomet- 
rical effect, and is also observed in lower energy mea- 
surements pi 124025!] . The multiplicity dependence of Ri 



differs from observations at lower energies; E895 |l4j and 
NA49 |2JJ observe no dependence of Ri over a wide range 
of multiplicity at the AGS and SPS, while NA44 (26) re- 
ports a sharp increase for the very highest multiplicity 
collisions at the SPS. 

The right panels of Fig. @ show, for the events in mul- 
tiplicity bin 3, the dependence of the HBT parameters on 
A increases with Tot, consistent with 
in which the increase 
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m T = \Jvt + m 

studies at lower energy 

was attributed to decreased contributions of pions from 

long-lived resonances at higher px- 

The radius parameters decrease significantly with my. 
This m-r-dependence, weak at ^/snn ~ 2 GcV [|14j and 



growing stronger with collision energy [£2-24 27]], reflects 
pion emission from a radially expanding source Jy. The 
mr-dependence at RHIC is similar, but not identical, to 
that observed in central Pb+Pb collisions at the CERN 
SPS |2^J2^,Q7|. The significant decrease of R with tut 
contrasts with the mT-independence of R measured at 
midrapidity in collisions at SPS energy 23 27|], For our 
highest rriT data, R < R s , possibly indicating significant 
source opacity resulting from the high particle densities 
generated in RHIC collisions p8| , p9[ . 

To place our results in context, Fig. shows the world's 
dataset of correlation parameters for midrapidity, low- 
p T (~ 170 MeV/c) tt~ from central Au+Au or Pb+Pb 
collisions. The parameter A falls smoothly and rapidly 
from unity (the expected value from naive assumptions) 
at y/sNN ~2 GeV, to about 0.5 at RHIC; this decrease is 
attributed partially to an increased fraction of 7r~ arising 
from long-lived decays as ^/snn increases H,p0[ • Since A 
is also affected by several experiment-specific effects (e.g. 
e~ rejection efficiency), we do not focus on the details of 
its excitation function. 

The radius parameter R s correlates most directly with 
source geometry Q|. After an initial decrease at AGS 
energies, attributed to increasing space-momentum cor- 
relations H] , R s appears to rise slightly with collision en- 
ergy, reflecting a larger freeze-out volume with increasing 
pion multiplicity. Similarly, Ri does not exhibit a large 
increase with collision energy, after the initial increase 
between AGS and SPS energies. 

The parameter R encodes both geometry and 
timescale information; in the absence of flow and opacity 
effects, the emission timescale (duration of freeze-out) 
is given by r = \J R 2 — R%//3t, where (3t is the av- 
erage transverse velocity of a pion pair. Both R and 
\JR 2 — R^ are relatively independent of energy; we see 
no evidence for a large increase in emission timescale (e.g. 
due to QGP formation ||). 

In hydrodynamic models flqfl, the ratio R /R s is a 
sensitive probe of the Equation of State (EOS). In a 
purely hadronic (non-QGP) scenario, R /R s « 1.0 — 1.2, 
while formation of QGP is predicted to produce R /R s « 
1.5 — 10. This ratio, plotted in the bottom panel of 
Fig. [| for low-pr pions, does not show a significant rise 



at RHIC. Furthermore, at RHIC R /R s is observed to 
decrease (from 1.07 to 0.89) with pt, in contrast to re- 
cent transport calculations pW which include effects of 
hadronic rescattering after a hydrodynamic stage with a 
QGP phase. We note, however, that the long emission 
duration QGP signature is expected || to vanish again 
at energies higher than its onset, due to rapid expansion 
of the system. Because of the large energy difference 
between SPS and RHIC, without intervening measure- 
ments, we cannot rule out an increase in the emission 
timescale at some lower, intermediate energy. 



In conclusion, STAR has measured two-pion correla- 
tion functions in Au+Au collisions at ^snn = 130 GeV. 
Transverse HBT radii grow with event multiplicity, re- 
flecting the evolution of source geometry with centrality. 
In contrast to studies at lower energies, Ri also increases 
steadily over a large range of multiplicities; more study is 
required to understand the origin of this effect. The m-r- 
dependence of the HBT radii is stronger than at lower en- 
ergies, suggesting emission from a more rapidly expand- 
ing source at RHIC, and in particular the much stronger 
decrease in R at high tot may indicate the onset of 
opaqueness in the dense system. The pion interferometry 
excitation function for the heaviest ions now spans nearly 
two decades in ,/snn- No sudden jumps in HBT radii 
are observed, but lower energy RHIC measurements are 
needed to complete the search for a predicted increase in 
emission timescale related to the possible onset of QGP 
formation. 
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FIG. 1. Coulomb-corrected (full dots) and uncorrected 
(circles) correlation functions for 1ow-pt 7r~ emitted at midra- 
pidity from central collisions. Shown in panel (a) is the Qi nv 
correlation function, and in panels (b-d) projections of the 
3-dimensional correlation function onto the q , q s , and qi axes. 
To project onto one (/-component, the others are integrated 
over the range 0-35 MeV/c. Fits to Coulomb-corrected cor- 
relations are shown by curves. 
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FIG. 2. The multiplicity dependence of the HBT param- 
eters is shown for low-pr pions on the left; bin 3 contains 
the high-multiplicity events. For high-multiplicity collisions, 
the mr-dependences are shown on the right. Error bars in- 
dicate statistical uncertainties; systematic uncertainties are 
represented by the height of the shaded regions below. 
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FIG. 3. The energy dependence of -k HBT parameters 
for central Au+Au (Pb+Pb) collisions at midrapidity and 
p T « 0.17 GeV/c Q||,||g||j. The SPS data are offset 
slightly in ^/snn for clarity. Error bars on NA44, NA49, and 
STAR results include systematic uncertainties; error bars on 
other results are statistical. 
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